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ABSTRACT 

Many proto-planetary nebulae (PPN) appear as narrow coUimated structures 
sometimes showing multiple, roughly aligned lobes. In addition, many PPN 
flows have been shown to have short acceleration times. In this paper we explore 
whether jet or "bullet" (a massive clump) models fit the observations of individual 
collimated lobes adequately by comparing simulations of both radiatively cooled 
(stable) jets and bullets. We find that the clump model is somewhat favored 
over jets because (1) it leads to greater collimation of outfiows (2) it accounts 
better and more naturally for ring-like structures observed in the PPN CRL 
618, and (3) it is more successful in reproducing the Hubble-fiow character of 
observed kinematics in some PPN. In addition, bullets naturally account for 
observed multipolar fiows, since the likely MHD launch mechanisms required to 
drive outfiows make multiple non-aligned jets unlikely. Thus we argue that PPN 
outfiows may be driven by explosive MHD launch mechanisms such as those 
discussed in the context of supernovae (SNe) and gamma-ray bursts(GRB). 

Subject headings: ISM:jets and outfiows — planetary nebulae:general — planetary 
nebulae individual (CRL-618) — stars:AGB and Post-AGB 
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1. Introduction 

In the past decade, images of very young planetary nebulae (PNs) and proto-planetary 
nebulae (PPNs) have revealed an unexpected diversity of morphological classes. Many of 
these objects appear to exhibit a level of complexity that cannot be accounted for in terms 
of the Generalized Interacting Stellar Winds model (GISW; Balick & Prank 2002 and ref- 
erences therein). Of particular interest are objects exhibiting point-symmetric, multi-polar, 
and "butterfly" morphologies, as well as bipolar and multi-polar objects exhibiting highly 
coUimated "jet-hke" outflows. 

The appearance of these coUimated and sometimes multi-polar outflows in so many 
PPNs has led to the suggestion that high-speed jets operate during the late asymptotic giant 
branch (AGB) and/or post-AGB evolutionary phases of the central star (Sahai Sz Trauger 
1998). While the GISW model can account for narrow jets (Icke et al. 1992; Mellema & Frank 
1997; Borkowski, Blondin, & Harrington 1997), it assumes the winds are radiatively driven. 
Radiative acceleration cannot however account for these flows since a number of observational 
studies demonstrate a momentum excess such that a factor ~ 10^ exists between outflow 
momentum observed and what can be attributed to stellar radiation pressure (Bujarrabal 
et al. 2001). Moreover, it is difficult to attribute the degree of observed coUimation to a 
large-scale dust torus as is usually required in the GISW model. In addition, the problem of 
accounting for the precession necessary for the production of point-symmetric flows remains 
to be solved [for an accretion disk based model see Icke (2003)]. For these reasons, the 
suggestion has been made that PPN jets and coUimated outflows are magnetically driven 
(Blackman et al. 2001a, 2001b; Frank & Blackman 2004; Matt, Frank, & Blackman 2006, 
Frank 2006). Magnetically driven models couple rotation to a magnetic fleld. Jets therefore 
are bound to flow along the rotational axis of the central object and it is difficult to see how 
multiple jets of similar size can be driven by such a mechanism. We discuss these models 
and this issue in more detail at the end of the paper. 

Observationally, clumps and coUimated flows occur in many stellar outflows though not 
always together. The outflows in Wolf-Rayet (WR) nebulae are clumpy, but jets are not 
observed. In young stellar objects (YSOs), jets and coUimated bipolar outflows are quite 
common, and while they can often be clumpy, the jet beams — distinct from the bow shocks 
which they drive — are often apparent, stretching all the way back to the stellar source. In 
mature PN, clumps are often seen [as in NGC 2392 (Eskimo), 6853 (Dumbbell) and NGC 
7293 (Helix)]. Fully articulated jets are however very rare. We note that ionization shadows 
and "mass loaded" flows behind clumps can give the appearance of jets. In some mature PN 
such as the Cats Eye nebula, structures appear (some of which fall under the term FLIERS) 
which may be the remnants of poleward-directed flows. In HST images of many PPNs, the 
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outlines of reflected light are often bipolar, but within these boundaries the illuminated gas 
seems irregular. Thin jets (as opposed to thin finger-shaped lobes) are rarely seen directly 
except (perhaps) in 0H231.8+04.2 (Calabash) . However, pairs or sets of knots lying along 
or near the apparent symmetry axes are not unusual (Ml-92, IRAS 20028 +3910, IRAS 
16594-3656, Hen 3-1475). Thus the creation of continious jets as in the case of YSOs does 
not seem to be the norm in PN and PPN. 

Clumps or "bullets" driven into the surrounding media have been found to be an effective 
explanation for some stellar outflow structures. In Poludnenko, Frank, & Mitran (2004) the 
authors modeled the strings of r] Car as bullets of high speed material ejected by the star. 
The simulations showed that long, thin morphologies similar to jets were readily obtained 
along with multiple rings associated with vortex shedding and the break-up of the clump. 
The authors suggested that such "impulsive" models may be useful in PPNs as well. Such 
a scenario is very different from the jet-driven explanation for PPN/PN. In this paper we 
seek to explore the usefulness of the clump picture. 

Soker (2000) has analytically explored the role of jets in PNs. In the excellent study 
of Lee & Sahai (2003), simulations of jets as the drivers of PPNs were presented, including 
detailed comparisons with observations. Our goals in this paper are more modest. In what 
follows we take a first step in the exploration of the clumps vs. jets issue by examining 
2-dimensional and 3-dimensional pairs of simulations, with each pair consisting of either a 
steady jet impinging upon a circumstellar gaseous medium, or a clump of gas which is fired 
ballistically through the same medium along a trajectory corresponding to the direction of 
flow in the steady jet 0. Both the jet and the clump are assumed to be magnetically launched 
though no attempt to model the launch mechanism is made here, and the simulations are 
purely hydrodynamic. For the present, we are merely interested in examining how the clump 
and jet differ in their effect on the surrounding circumstellar medium. As we will show, the 
jet and clump models show differences which require further study, but the clumps provide 
at least as good, or better, an account of key observational characteristics. Given the fact 
that in some cases multiple outflows are seen in a single object (such as CRL 618), the clump 
model may be more plausible since what are often interpreted as multiple "jets" could instead 
arise naturally from the fragmentation of an explosively driven polar directed shell. We note 
here that none of the widely accepted magnetically launched outflow models would create 
continuous multiple jets of similar or equal age driven in slightly different directions. 

We note also that new models of binary stars in the context of PN's (Nordhaus & 



^In a related study Raga et al.(2007) have also recently presented a model of the "3D structure of a 
radiative, cosmic bullet flow." 
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Blackmail 2006, Nordhaus, Blackmail & Frank 2007) show the extent to which envelope 
ejection can be shaped by gravitational interactions. In Nordhaus & Blackman 2006 Common 
Envelope scenarios which lead to aspherical mass loss (including disk creation and possible 
MHD launching) were articulated. In Nordhaus, Blackman & Frank 2007 Common Envelope 
models were explored as the source of differential rotation in the primary which could drive 
strong dynamo supported magnetic fields. These models showed that while single stars may, 
in some cases, be able to support a strong field over AGB timescales, binary interactions 
were highly effective at creating the fields needed to power PPN outflows at the evolutionary 
moment when they will be required. As we will see, such models provide strong theoretical 
support for the scenario we argue for in this paper. 

In section |2] we provide information concerning the numerical methods used, details of 
the jet and clump models, and a discussion of the initial conditions. In section [3] we discuss 
the results of our simulations and in section H] we summarize our conclusions. 



We have carried out two pairs of hydrodynamic simulations (one "medium-resolution" 
3D pair and one "high-resolution" 2.5D pair). Each pair consists of a jet and clump re- 
spectively with each parameterized to be as similar to one another as possible. Specific 
parameter values for the jet, clump, and ambient medium for each simulation are given in 
table [H The simulations are performed using the AstroBEAR code which is an extension of 
the BEARCLAW adaptive mesh refinement (AMR) package for solving conservation laws. 
[For a detailed description of the AstroBEAR package see section 3 of Cunningham Frank 
& Blackman (2006).] The domain is a rectangular box with a square cross-section and with 
the X-axis chosen to intersect the center of the left square face of the domain. The origin 
of coordinates is placed at this point of intersection. The clump and jet are launched along 
the X-axis and placed so that their centers coincide with it. The jet was modeled in 3D with 
a circular cross-section of radius rj (in 2.5D the jet cross-section reduces to a line-segment 
of length 2rj and a thickness of one computational cell). The jet was launched into the 
domain from a set of fixed cells along the domain boundary. To prevent the expansion of 
the jet inflow boundary with time, a ring of zero velocity and with outer radius 1.125ro was 
maintained around the jet-launching region. The velocity profile of the jet was smoothed 
about a nominal value Vj^ according to 



2. Computational Methods and Initial Conditions 




(1) 
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where s is a shearing parameter taking values between and 1 and for the jet simulation 
presented here is set equal to s = 0.9. The clump was modeled in 3D as a spherical over- 
density of radius = Vy (The sphere reduces to a circle in 2.5D.) Its initial position in the 
domain is chosen so that its center is located at the point 

re,o(x,z/,^) = (2re,0,0), (2) 

and the density of the clump as a function of location within the clump is 



nc{r) = n„(r) + Uq 

where 



2 

\r - rc,o 

To 



(3) 



n^{r) = min I hq, \ , (4) 

is the ambient number density profile in regions of the domain unoccupied by jet or clump 
gas, and where — x'^ + + , no is the nominal ambient number density, and ro is a 
characteristic length taken to be equal to the jet or clump radius. The 3D (2.5D) simulations 
are carried out on a base grid with a resolution of 6 (12) cells per jet/clump radius and with 
two levels of AMR refinement providing an effective resolution of 24 (48) cells per jet/clump 
radius. In all cases, radiative cooling is modeled using the atomic line cooling function 
of Delgarno & McCray (1972), and we do not attempt to follow the detailed ionization 
dynamics or chemistry of the cooling gas. Given that both models give rise to similarly 
expanding shells of shock-heated gas we do not expect this simplification to materially affect 
our conclusions. 
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Model 


Pcirameter 


Value (2.5D) 


Value (3D) 


Jet 


Radius, rj 


500 AU 


500 AU 




Computational cells per rj . 


48 


24 




number density, Uj 


500 cm-3 


500 cm-3 




peak velocity, Vj q 


100 km s^l 


100 km s-l 




Temperature, Tj 


200 K 


200 K 




Nominal ambient density, Ua 


500 cm^"^ 


500 cm-3 




Ambient temperature. To . . 


200 K 


200 K 




Shear parameter s 


0.9 


0.9 


Clump 


Radius, r'c 


500 AU 


500 AU 




Computational cells per Tc ■ 


48 


24 




nominal number density, rio 


500 cm-3 


500 cm-3 




velocity, Vc 


100 km s-l 


100km s-i 




Temperature, Tc 


200 K 


200 K 




Nominal ambient density, Ua 


500 cm-3 


500 cm-3 




Ambient temperature, Ta ■ . 


200 K 


200 K 



Table 1: Simulation Parameters 
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3. Results: Morphology 

Results of our simulations are presented in figures [Tl-|6l In Figures [1] and [2] we present the 
results of medium resolution, (24 cells per radius), 3D simulations of one jet and one clump 
respectively. The length of the domain in these simulations is 20 computational units with 
one computational unit corresponding to a physical scale of 500 AU. (One computational 
unit is also the value chosen for the radii of the jet and clump.) In each figure the upper 
image shows a plot of emission integrated along the line of sight which in the case of these 
figures is perpendicular to the plane of the image. The lower image in each figure is a plot 
of the logarithm of density in a plane coincident with the x-y plane. The images show the 
jet and clump near the end of their respective runs at time t ~ 498 yr for the case of the 
clump and at time t ~ 636 yr for the case of the jet. The resolution is seen to be sufficient 
to capture vortex-shedding events in both simulations. It is also evident from these images 
that while one can easily distinguish jet from clump in the density maps, the emission maps 
are quite similar. We note that the clump gives rise to a somewhat more coUimated flow, 
while the jet bow shock expands laterally at a greater rate than the clump bow shock. The 
jet also lags behind the clump in its forward motion. This is likely due to the streamlining 
that occurs as the head of the clump is reduced in size as material is ablated away via its 
interaction with the ambient medium. 

Due to limits on computational resources, it was necessary to impose limits on the 
resolution and run time of the 3D simulations from which the images in flgures [1] and [2] are 
taken. The simulations end just as the vortex-shedding events begin to have an interesting 
effect on the nebular environment. To explore this stage further we carried out the second 
pair of 2.5D simulations mentioned above. In these high-resolution simulations, the effective 
resolution was doubled to 48 cells per jet/clump radius, the length of the domain was doubled, 
and the transverse dimensions of the domain were enlarged in an attempt to accomodate 
the lateral expansion of the jet/clump bow shock (this latter adjustment was successful only 
for the case of the clump). Results from these simulations are presented in flgures [3] through 
El In flgures [3] and H] we again present images of the logarithm of density for the jet and 
clump respectively — this time reflected about the axis of symmetry. Both flgures show the 
simulations at various stages of the flow. We observe that the differences found between the 
two models in the 3D simulations — i.e. the faster domain crossing time and the higher degree 
of coUimation exhibited by the clump — are seen again in these images. The vortex shedding 
however, is now captured with greater clarity for both jet and clump, and we begin to see 
signiflcant qualitative differences in the manner in which these events unfold. In particular 
we note that shedding events are much more frequent in the case of the clump. These results 
mirror those found by Poludnenko, Frank & Mitran (2004). It is noteworthy that their study 
used a different integration scheme than used here. AstroBEAR has a number of schemes 
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built into it and in the Poludnenko study a Wave Propagation scheme was used (LeVeque 
1997), while here a MUSCL-Hancock method is used. The fact that the basic morphology 
of clumps driving bow shocks dominated by vortex shedding events is recovered using both 
schemes gives us confidence in this aspect of the dynamics. 

3.1. Morphology 

To get a better sense of how the differences between the models might appear observa- 
tionally, we present integrated emission maps for the jet and clump respectively in figures 
E] and O These figures were produced by calculating the effect on the line-of-sight emission 
resulting from a rotation of the cylindrically symmetric data set about the axis of symme- 
try. The intensity shown, which does not distinguish among cooling lines, was determined 
according to: 

Ii,j,k = ^knlj^kAiTij^k), (5) 

where and k refer to the x,y and z directions in the final data cube created by rotating 
n{r, z) and T(r, z) about the axis of symmetry, and A is the cooling function. The images 
shown correspond to the final frames in each of figures [3] and H] respectively. Each figure 
provides two views of the data: one in which the angle of inclination of the symmetry axis 
with respect to the image plane, 6', is 0°; and one in which it is 20°. One difference between 
the jet and clump cases appears in the shape of the head of the bow shock. A clump has a 
finite reservoir of mass which interacts with the ambient medium. As the clump propagates 
down the grid, it drives a (bow) shock wave into the ambient medium. A second shock passes 
through the clump heating and compressing it. When cooling is present this "transmitted 
shock" first leaves the clump flattened. As material is then ablated away via the interactions 
with the ambient medium the remaining clump material becomes dense and streamlined in 
the direction of propagation. At later times in the simulation the dense core of the clump 
drives a V - shaped bow shock head. In the case of a jet the situation is different. The jet 
head drives a bow shock into the ambient medium and transmitted shock, called a jet shock, 
propagates back into the jet material. Decelerated jet material flows transverse the these 
shocks inflating a cocoon behind the wings of the bow shock. Unlike the clump however, 
there is always more high speed material behind the jet shock/bow shock pair to resupply 
the interaction. Thus with material continuously flowing into the cocoon, the bow shock 
head remains wider and takes on a flatter more U - shaped configuration. Such a distinction 
between V - and U - shaped flows may be important in comparing with observations. We 
note that both the 2.5D and 3D simulations both show this difference. We note however 
that the axial symmetry will tend to enhance features on the axis. Our 3D runs do not yet 
have the resolution to accurately track the break-up of the clump. Thus the V and U bow 
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shock head distinction must be considered less than conclusive and await further study. 

Vortex shedding provides another morphological distinction. In the case of the clump, 
the relatively frequent shedding events have led to a series of thin, irregularly spaced rings of 
enhanced intensity centered about the symmetry axis. These are reminiscent of the ring-like 
structures observed in some coUimated PPN outflows [see for example Trammell & Goodrich 
(2002)]. The shedding events occurring in the jet simulation lead to similar structures, but 
these are less frequent and somewhat more band-like in character. The qualitative differences 
in the manner in which these rings form in the outflows depending on whether one models 
them as jets or clumps might suggest a means of distinguishing between the two models in 
observations. One must, once again, be careful not to over- interpret these results due to 
limits on the resolution and the fact that these simulations are 2.5D. We thus conclude that 
the high-resolution 2.5D simulations lend weight to the assertion that clumps and/or jets 
can account for observed ringed structures, but neither can be ruled out as a model for the 
collimated outflows observed in the environments of PPNs. In the meantime, we note that 
this conclusion in itself is important with respect to PPN studies as we will discuss in the 
last section. 

It is also noteworthy that Lee 8z Sahai (2003) attempted to model the rings via a pulsed 
jet. Each ring became associated with an "internal working surface" where faster moving 
material swept over slower moving material. The internal shocks lead to transverse motions 
of shocked material which impinge upon the bow shock. As might be expected, the strength 
of the emission from these shocks decreased as the pulse traveled down the length of the 
beam. Such dimming of the rings with distance from the source is not what is observed in 
CRL 618. The clump on the other hand produces the opposite kind of pattern, as ablation 
events on the clump lead to rings that are bright closer to the head of the bow shock. 

3.2. Kinematics 

In addition to comparisons of morphology, it is also important to consider the flow 
kinematics, since observed flows are often seen to exhibit "Hubble-flow" characteristics — 
that is, the flow velocity is observed to increase linearly with respect to distance from the 
flow origin. To address this issue we present in figure [TJ plots of the x-component of the 
flow velocity averaged over the directions transverse to the flow (vx). Velocity tracers were 
not used in our simulations. In order, therefore, to differentiate between mildly perturbed 
ambient gas, and gas that is fully involved in the flow, values of velocity < O.Olfo were 
ignored, where vq is the initial velocity of the clump or jet gas. The top row of the figure 
shows, from left to right, the results for the 3D jet and clump respectively while the bottom 
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figure shows the results for the 2.5D jet and clump. The data for these plots are taken from 
times near the end of the simulation when the flows have crossed most of the domain. In 
the two plots involving clumps, there are large regions of the flow for which the variation 
of velocity with distance is roughly linear. The jets on the other hand, fail to model this 
behavior altogether. Comparison of the 3D clump plots to the 2.5D case yields an interesting 
result. Recall that the first vortex-shedding events in both clump and jet were observed to 
occur in the 3D simulations shortly before the end of the simulation. Because of this they do 
not have time to perturb the flow in a way that might be noticeable in these plots. However, 
when we examine this phenomenon kinematically in the extended spatial domain allowed by 
the 2.5D simulation, we find that while the vortex shedding events, appear to perturb the 
kinematics of both jet and clump, these perturbations do not alter the overall qualitative 
character of the flows in either case. 

To further examine the kinematics of our simulated flows, we have also produced a set 
of synthetic position-velocity (PV) diagrams for the 2.5D simulations. These are presented 
in figure [HI Once again, as in figures [5] and [6] we present our results in pairs corresponding to 
values of 0° and 20° for the angle of inclination, 6', of the flow symmetry axis with respect to 
the image plane. These plots were produced by calculating the velocity structure along the 
line of sight and with the "slit" placement taken to be along the projected axis of symmetry. 
Results for the jet are given in the first row of the figure, and results of the clump are given 
in the second row. In these images, the difference between the clump and jet are even more 
striking. For either angle of inclination, the velocity structure of the jet cannot be said to 
be even approximately Hubble-like. The clump however, continues to exhibit line-of-sight 
velocity structure indicative of a linear increase with distance along the projected direction 
of flow. The effect is particularly apparent in the case of the flow which is inclined with 
respect to the image plane. These results, and those of figure [7] suggest that it may be 
possible to distinguish outflows from steady jets from those from explosive events through a 
careful examination of their kinematics. 

One caveat which must be considered in these results is the role of emission. Because 
our models do not track emission from individual species we cannot separate the emission at 
the bow shock from that within the jet or from the shocked jet material. In Ostriker et al. 
(2001), a model for the emission from a jet-driven bow shock was presented which showed 
a characteristic spur pattern in synthetic PV diagrams. The spur exhibits a rapid drop in 
velocity away from the tip of the bow shock. Lee & Sahai (2003) found a range of patterns 
in their jet simulations which in some cases took on the spur morphology. Thus our results 
are suggestive of the differences between jets and clumps and indicate that clumps appear 
to be better, in general, at recovering quasi-linear increases in velocity along the nebular 
outflow lobe. 



3.3. Kinematic Models 



In order to interpret our results we consider the time-dependent distortion of the clump 
gas during the evolution of the outflow. Strongly radiative, hypersonic clouds of any geometry 
will be rapidly compressed into a thin ballistic sheet after ejection by the outflow progenitor. 
We therefore consider the motion of a cylindrically symmetric disk with surface density x(r) 
and velocity v{r,t) where v{r,0) — Vo to model the time-dependent evolution of the clump 
gas. The equation of motion for a differential ring of the disk under the ram pressure of the 
ambient gas of density p is given by: 

pv {r,t) = -X-^- (6) 

Because the outflow bow shock is convex, most of the outflow-entrained ambient gas will be 
swept outside of the path of the clump into the bow shock. We therefore neglect accretion 
of ambient material onto the clump and the kinematics of ambient material ejection in this 
model. Because the disk is hypersonic in a strongly cooling environment, we consider the 
model disk to be ballistic, neglecting pressure forces. For simplicity we also take the density 
of the ambient gas to be constant. Thus the equation of motion for a differential ring of 
clump gas with radius r integrates to: 

V — — . (7 

l + pvot/x 

The distance traversed by the ring is given as: 



L{x,t)^ [ v{t)dt^^ In 
Jo P 



1 + 

X 



(8) 



The quantities v, t, and L, all refer to the same ring. What differentiates one ring of 
material from another is the parameter r. Now at some late time t, we imagine the rings 
to have been distributed over the length of the outflow with this distribution depending on 
r. For this flxed value of t, we are interested plotting the velocity of each ring against its 
corresponding distance. The r-dependence of v and L enter into the expressions for these 
quantities through the surface density x- We therefore model x by assuming that the clump 
of gas from which our disk formed was initially spherical, of constant volume mass density 
a, and compressed in such a way that all material within the volume of the clump and lying 
along a given line passing through the clump in the direction of its motion remains on this 
line after compression. Then, 



X(r) = 2aro\/T^^Jr/roj^, 



(9) 
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where tq is the radius the clump/disk. Introducing the dimensionless quantities: 



r 

r = — , 

V 

V = — 
r pVot 



2(Tr, 







and, 



our parametric equations are: 



2aro 



1 + 



and 



L = Vl - r2 In 



1 



(10) 

(11) 

(12) 
(13) 

(14) 
(15) 



The value of t is chosen by assuming cr > 2p, and by noting that at late times VqI > tq. 
In our 2.5D simulations we have vqI/tq ~ 4t ~ 20 making t ~ 5. Figure [9] shows a plot of 
f vs L in astronomically relevant units with this choice of t. For purposes of comparison, 
a line with an appropriately chosen slope and intercept is plotted as well. In spite of the 
simplicity of our analytical model, a comparison of this plot with the upper-right-hand plot 
of figure [7| reveals good agreement between the two. Both curves exhibit a small concave 
curvature for small L while becoming increasingly linear with increasing L (the downward 
turn in the 2.5D simulation-based plot results from extending the plot into regions not yet 
reached by the clump). We also find that the range of L- values over which the curve shows 
linear behavior increases for increasing t, implying that a correlation between the kinematical 
ages of PPN outflows and the extent to which they are observed to exhibit Hubble-fiow. We 
take this calculation as further evidence that our simulations are accurately capturing the 
dynamics of the flow and the relevance of bullet models for PPN. 
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Fig. 1. — 3D Jet at time t ~ 636 yr, top: integrated emission assuming atomic line cooling; 
bottom, base-ten logarithm of density. 
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Fig. 2. — 3D clump at time t ~ 498 yr, top: integrated emission assuming atomic line 
cooling; bottom, base-ten logarithm of density. 
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Fig. 3. — Density of 2.5D Jet at times t ~ 336, 628, 896, and 1165 yr. The data is shown 
here reflected about the axis of symmetry. 
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Fig. 4.— Density crosscut of 2.5D Clump at times t ~ 208, 477, 753, and 1082 yr. The data 
is shown here reflected about the axis of symmetry. 
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Fig. 5. — Integrated emission of 2.5D Jet at time t ~ 1132 yr. The data is shown here 
rotated about the axis of symmetry and with angles of inchnation of the symmetry axis with 
respect to the image plain of ^ = O°(top), and 9 — 20°(bottom). 
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Fig. 6. — Integrated emission of 2.5D Clump at time t ~ 1082 yr. The data is shown here 
rotated about the axis of symmetry and with angles of inclination of the symmetry axis with 
respect to the image plain of ^ = O°(top), and 9 — 20° (bottom). 
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Fig. 7. — Comparisons of weighted average parallel flow velocity, {vx), vs distance from the 
flow origin for clump and jet for the 3D simulations (top row) calculated at time t ~ 636 yr, 
and the 2.5D simulations(bottom row) calculated at time t ~ 1165 yr. (See the text for an 
explanation of the weighting.) 
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Fig. 8. — Position-velocity diagrams at time t ~ 1082 yr for the 2.5D Jet (top row) and the 
2.5D Clump (bottom row) assuming indination angles of ^ = 0°(left) and 9 — 20° (right). 
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Fig. 9. — V{r) vs L{r) for t = 5. See section [331 for an explanation of this figure. 
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4. Discussion and Conclusions 

In this paper we have examined the resuhs of two pairs of simulations intended to model 
the gross morphological and kinematical properties of PPNs. Our primary purpose was to 
ask not if we could distinguish between jet and clump models, but instead to ascertain if 
clump models could perform equally well at recovering these properties. Below we explain 
the justification for this more explicitly and present the conclusion of our study. As was 
discussed in the introduction, MHD models of PN shaping have been explored by a variety 
of authors and in a variety of forms. As we will explain below, it is the imperative of MHD 
models, particularly magneto-centrifugal launch and coUimation scenarios, which motivate 
this paper. 

Two distinct classes of model for the magnetic shaping of winds in PNs and PPNs have 
been suggested to date. First there is the Magnetized Wind Bubble (MWB) model, originally 
proposed by Chavalier & Luo (1994) and studied numerically by Rozyczka & Franco (1996) 
and Garcia-Segura et al. (1999). In these models, an initially weak toroidal magnetic field 
is embedded in a radiatively-driven wind. This configuration has been shown capable of 
accounting for a wide variety of outfiow morphologies including highly - collimated jets. 
These models clearly demonstrate the importance of magnetic fields. They cannot however 
account for the excess momentum in the flows (Bujarrabal et al. 2001) because of the weak 
flelds which simply ride along in a radiatively-driven wind. 

The other class of model invokes so-called Magnetocentrifugal Launching (MCL; Bland- 
ford & Payne 1982; Pelletier & Pudritz 1992). This paradigm has, for many years, been 
explored as the mechanism driving jets in young stellar objects (YSOs), micro-quasars and 
active galactic nuclei ( AGNs) . The MGL paradigm assumes the presence of a rotating central 
gravitating object (which may or may not include an accretion disk). In the case of a disk, 
plasma is threaded by a magnetic field whose poloidal component is in co-rotation with the 
disk. Disk-coronal gas is then subject to centrifugal force which accelerates the gas hinging 
it out along field lines. The magnetized plasma eventually expands to a configuration where 
the toroidal component of the field dominates and hoop stresses collimate the fiow. Thus 
the MGL paradigm accounts for both the origin of the wind and the means of coUimation. 

The success of the MGL paradigm in modeling jets associated with YSOs and AGNs has 
led some authors to suggest applying the idea in the context of PNs and PPNs (Blackman et 
al. 2001a, Frank & Blackman 2004,). Most recently it has been shown that the observed total 
energy and momentum in PPNs can be recovered with disk wind models using existing disk 
formation scenarios via binary interaction (Prank & Blackman 2004 and references within). 

Most theoretical investigations of the MGL paradigm assume a steady-state fiow. Ob- 
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servations suggest however, that acceleration times for the flows are as much as an order of 
magnitude shorter than typical kinematical PPN ages (Bujarrabal et al. 2001). This implies 
that the mechanism responsible for the observed flows may operate explosively, i.e. the time 
over which the mechanism acts is short compared to the lifetime of the flow. Moreover, it 
has been suggested by Alcolea et al. (2001) that such a scenario would also provide the 
most straightforward explanation for the "Hubble law" kinematics observed in some PPN 
outflows (Bahck & Prank 2002; Bujarrabal, Alcolea, & Neri 1998; Olafsson & Nyman 1999). 
The MCL paradigm can act transiently however, when linked with the rapid evolution of its 
source, as for example in the case of the proposed mechanisms for gamma-ray bursts (GRBs; 
Piran 2005) and Supernovae (SNe). This scenario has been investigated by a number of 
authors (Kluzniak & Ruderman 1998; Wheeler et al. 2002; Akiyama et al. 2003; Blackman 
et al. 2006). In these scenarios differential rotation twists an initially weak poloidal flcld 
thereby generating and amplifying a toroidal flcld. When the toroidal component reaches a 
critical value it drives through the stratifled layers of the collapsing core carrying trapped 
material with it. The hoop stresses associated with such a fleld also serve to coUimatc the 
flow. 

Recently Matt, Frank, & Blackman (2006) examined numerically a simplifled version of 
this idea which was originally suggested for PPN in Blackman et al. 2001b. In these studies 
the authors began with a gravitating core threaded by an initially poloidal fleld set rotating 
at 10% of the escape speed within an envelope of ionized gas. As the simulation progressed, 
the resulting toroidal fleld was sufficiently strong to drive a complete and rapid expulsion 
of the gaseous envelope. Since the initial conditions assumed in Matt, Prank & Blackman 
(2006) are applicable to either a young PPN or a collapsing protopulsar, it is reasonable to 
ask whether such transient events are occurring in the early stages of the formation of PNs, 
and if such events can serve as well as steady-state jets can in accounting for the complex 
morphologies observed in such systems. We note that these classes of model arc sometimes 
referred to as "magnetic towers" or "springs" because it is the gradient of toroidal fleld 
pressure which drives the outflow. Again we note that the magnetic flelds needed for our 
scenario can be delivered by binary interactions as has been demonstrated by Nordhaus & 
Blackman 2006 and Nordhaus, Blackman & Prank 2007. 

There is growing evidence to suggest that magneto-centrifugal launch models are appro- 
priate for PNs and, more importantly, PPNs (Vlemmings et al. 2006). Taken together with 
the evidence that many PPNs have short acceleration time scales for which Tqcc < -iTdyn it 
suggests that some PPNs may be considered to have arisen from explosive, or at least im- 
pulsive, depositions of momentum and energy into surrounding, circumstellar environments. 
Together we argue that these lines of evidence suggest that many PPNs may be shaped by 
shells which fragment into clumps rather than multiple jets. 
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There is a subset of PPNs and young PNs with multiple lobes of roughly similar size. 
These include CRL 2688, CRL 618, IRS 19024+004, IRS 09371+1212, Ml-37 and He2-47. 
It is natural to try and interpret these structures initially as resulting from the action of jets. 
However, consideration of magneto-centrifugal launching models shows this to be unlikely. 
In all forms of the model, gravitational binding energy is tapped via rotational motions of 
the central source about some axis a;, and is converted into outflow kinetic energy using 
the magnetic fields as a "drive belt". The existence of a quasi-stable rotational axis is a 
requirement of the models in order to produce a continuous outfiow. Multiple jets of equal 
length are difficult to imagine in such a scenario as the jets would then each require their 
own rotational engine with separate alignments. Even so-called magnetic tower models which 
drive the jet by winding up an initially weak poloidal field require a net spin axis such that 
2Tin^Bp where is the number of turns about uj. 

The production of multiple bow shocks from clumps or bullets driven by a transient 
MCL process is not as difficult to envision. In Matt, Prank, & Blackman (2006), it was 
shown that the static envelope or atmosphere of a star could be entirely driven off of a 
rotating magnetized core. These models relied on the magnetic tower "spring" mechanisms, 
and the envelope becomes compressed into a thin shell which rides at the front of the ex- 
panding magnetic tower. Such a thin accelerating shell would be subject to a variety of 
instabilities including the Rayleigh- Taylor, Thin-Shell (Vishniac 1993) and Non-linear Thin 
Shell (Vishniac 1994) modes, all of which would be modified by the presence of an ordered 
magnetic field which would impose a long coherence length onto the resulting fiow. The 
precise detail of such fragmentation in this situation have yet to be calculated and stand as 
an open problem. Given the impulsive acceleration of a dense, radiatively cooling shell into 
a lower density environment, it is likely that the shell would fragment into a number of high 
mach number clumps directed along the poles (and perhaps the equator, see Matt, Frank, & 
Blackman 2006 and CRL 6888). A potential challenge to this model would be the creation of 
fragmentation modes that can, in some cases produce roughly equivalent clumps in terms of 
propagation direction on either side of the source as is seen in some cases. Given that caveat 
however, the ability for explosive magnetic tower models already explored in the literature 
to drive unstable shells makes them an attractive means of producing high mach number 
clumps. As we have shown, these clumps, propelled into the surrounding media then drive 
bow shocks which do at least as good a job as, if not better than, jets in recovering gross 
morphological and kinematic observations. Similar conclusions have recently been reached 
by Raga et al. (2007). 

In summary, the results presented here add weight to an emerging paradigm in which 
transient (explosive) MCL processes act as the driver for PPN evolution in some case. The 
fact that such magnetic launch mechanisms are already favored by some theorists to explain 
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supernovae and gamma-ray bursts (Piran 2005), makes all the more compelling the notion 
that lower energy analogues of the processes believed to be occurring during the penultimate 
stages of massive stars' evolution, are also occurring in low and intermediate mass stars. 
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